Single- qubit- gate error below 10 in a trapped ion 
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With a ^Be+ trapped-ion hyperfine-states qubit, we demonstrate an error probability per randomized single- 
qubit gate of 2.0(2) x 10^^, below the threshold estimate of 10^"* commonly considered sufficient for fault- 
tolerant quantum computing. The ^Be+ ion is trapped above a microfabricated surface-electrode ion trap and is 
manipulated with microwaves applied to a trap electrode. The achievement of low single-qubit-gate errors is an 
essential step toward the construction of a scalable quantum computer. 



In theory, quantum computers can solve certain problems 
much more efficiently than classical computers fV\. This has 
motivated experimental efforts to construct and to verify de- 
vices that manipulate quantum bits (qubits) in a variety of 
physical systems |2|. The power of quantum computers de- 
pends on the ability to accurately control sensitive superpo- 
sition amplitudes by means of quantum gates, and errors in 
these gates are a chief obstacle to building quantum comput- 
ers [3). Small gate errors would enable fault-tolerant oper- 
ation through the use of quantum error correction protocols 
m. While the maximum tolerable error varies between cor- 
rection strategies, there is a consensus that lO^'* is an impor- 
tant threshold to breach fHISl. Single-qubit gates with errors 
slightly above this level have been achieved with nuclear spins 
in liquid-state nuclear-magnetic resonance experiments |6| 
and with neutral atoms confined in optical lattices [7j; here we 
demonstrate single-qubit error probabilities of 2.0(2) x 10^^, 
substantially below the threshold. Reaching fault-tolerance 
still requires reducing two-qubit-gate errors from the current 
state of the art (7 x 10"^ for laser-based |i8l and 0.24 for 
microwave-based gates |9|) to similar levels. 

To determine the average error per gate (EPG), we use 
the method of randomized benchmarking ifTOl . Compared to 
other methods for evaluating gate performance, such as quan- 
tum process tomography 111], randomized benchmarking of- 
fers the advantage that it efficiently and separately can deter- 
mine the EPG and the combined state-preparation and mea- 
surement errors. Because it involves long sequences of ran- 
dom gates, it is sensitive to errors occurring when gates are 
used in arbitrary computations. In randomized benchmark- 
ing, the qubit, initialized close to a pure quantum state, is sub- 
jected to predetermined sequences of randomly selected Clif- 
ford gates 1 12 1 for which, in the absence of errors, the mea- 
surement outcome is deterministic and efficiently predictable. 
Clifford gates include the basic unitary gates of most pro- 
posed fault-tolerant quantum computing architectures. To- 
gether with certain single-qubit states and measurements, they 
suffice for universal quantum computing llT2l IT3l . To estab- 
lish the EPG, the actual measurement and predicted outcome 
are compared for many random sequences of different lengths. 
Under assumptions presented in Ref. 1 10|, this yields an aver- 
age fidelity as a function of the number of gates that decreases 
exponentially to 1/2 and determines the EPG. Randomized 
benchmarking has been used to quantify single-qubit EPGs in 
a variety of systems as summarized in Table |l] 



To improve on the results of Ref. ifTOl . we integrated a mi- 
crowave antenna into a surface-electrode trap structure |16|. 
The use of microwave radiation instead of optical stimulated- 
Raman transitions to drive qubit rotations suppresses deco- 
herence from laser beam pointing instability and power fluc- 
tuations and eliminates decoherence from spontaneous emis- 
sion. The microwave amplitude can be stabilized more eas- 
ily than laser power, and because the antenna is integrated 
into the trap electrodes, unwanted motion of the trap does 
not affect the microwave-ion-coupling strength. The small 
distance (40 jim) between the trap surface and the ion per- 
mits transition rates comparable to those based on lasers. 
Improved shielding from ambient magnetic-field fluctuations 
was achieved by locating the trap inside a copper vacuum en- 
closure held at 4.2 K by a helium-bath cryostat. The thickness 
of the walls, combined with the increase in electrical conduc- 
tivity of copper at 4.2 K, effectively shields against the ambi- 
ent magnetic field fluctuations that typically limit coherence 
in room-temperature ion-trap experiments ifTOl . This shield- 
ing is evident when we change the magnetic field external to 
the cryostat; the accompanying response in ion fluorescence 
lags the change with an exponential time constant of 3.8(2) s. 
In addition, cryogenic operation decreases the background gas 
pressure to negligible levels, thereby enabling long experi- 
mental runs with the same ion, and it suppresses ion heating 

EHm. 

The ^Be+ ion is trapped 40 jxm above a surface-electrode 
trap II20II constructed of 8-/im-thick gold electrodes electro- 
plated onto a crystalline quartz substrate and separated by 5- 
fim gaps (Fig. [TJ. A static magnetic field Bq, parallel to the 
trap surface and collinear with a Doppler cooling laser beam, 
is applied to break the degeneracy of the ground-state Zee- 
man sublevels (Fig. [T| inset). We drive 2s hyperfine 
transitions with microwave pulses near 1.25 GHz, coupled 
with a 4-nF capacitor to one end of a trap control electrode. 
The microwave current is shunted to ground at the other end 
of the electrode by the 4-nF capacitor of an RC filter. Mi- 
crowave pulses are created by frequency quadrupling the out- 
put of a direct-digital synthesizer whose frequency and phase 
can be updated in less than 1 /is by a field-programmable 
gate array [(FPGA), 16-ns timing resolution]. An rf switch 
creates approximately rectangular-shaped pulses. This sig- 
nal is amplified and is delivered via a coaxial cable within 
the cryostat and a feedthrough in the copper vacuum enclo- 
sure. In this Rapid Communication, we use the clock transi- 
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Reference 



System 



Gate error 



This Rapid Communication (2011) 
Reference [61 (2009) 
Reference (T) (2010) 
Reference (14| (2009) 
Reference fTOl (2008) 
Reference ti5J (2010) 



Single trapped ion 
Nuclear magnetic resonance 
Atoms in an optical lattice 
Trapped-ion crystal 
Single trapped ion 
Superconducting transmon 



2.0(2) X 10-^ 
1.3(1) X 10-4 
1.4(1) X lO"'* 

8(1) X lO-'* 
4.8(2) X 10"^ 

7(5) X 10-3 



TABLE I: Reported average EPG for Pauli-randomized n/2 gates in different systems as determined by randomized benchmarking. 
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FIG. 1: (Color online) Micrograph of the ion trap, showing radio- 
frequency (rf) electrodes and control electrodes. The red sphere indi- 
cates the approximate ion position in the x-y plane. Also shown are 
the directions of the static magnetic field Sq and of the microwave 
current used to drive hyperfine transitions. (Inset) Energy level di- 
agram (not to scale) of the 2s ^Si/2 hyperfine states in ^Be+. Blue 
dashed lines indicate the transitions used to prepare and measure 1 1) . 
The solid black line indicates the qubit transition, and the red dashed 
line indicates one of the transitions used to shelve | t) into a dark 
state. 

tion {\F = 2,OTf = 0) EE I ;) o |1,0) = I t)) in ''Be+ for the 
qubit instead of the previously used |2,— 2) — 1) tran- 

sition fTOl (Fig. [T] inset). The clock transition is a factor of 
20 less sensitive to magnetic-field fluctuations (950 MHz/T at 
field Bo = 1.51 X 10-3 T, compared to 21 GHz/T). 

A benchmarking experiment proceeds as follows. The ion 
is Doppler cooled and optically pumped to the |2,— 2) state 
with (7--polarized laser radiation near the pSi/2,2,— 2) f^- 
PP3/273,— 3) cycling transition at 313 nm. Then, the qubit is 
initialized in | J,) with two microwave % pulses, resonant with 
the |2,— 2) o |1, — 1) and |1,-1) O 1 1) transitions (blue lines 
in Fig. [T] inset). Pulse duration is then controlled by a digi- 
tal delay generator, which has a 5-ps timing resolution. The 
frequency, phase, and triggering of each pulse remain under 
control of the FPGA. 

A predetermined sequence of randomized computational 
gates is then applied. Each computational gate consists of 
a Pauli gate (n pulse) followed by a (non-Pauli) Clifford 
gate (7r/2 pulse). The gate sequence is followed by mea- 



surement randomization consisting of a random Pauli gate 
and a Clifford gate chosen deterministically to yield an ex- 
pected measurement outcome of either | f) or | J,) . The Pauli 
gates are chosen with equal probability from the set e^^^'^p/^^ 
where dp G {±0,, icT,., ±ov, ±/}. The Clifford gates are 
chosen with equal probability from the set where 
ffc £ ±<7v}- In practice, a Clifford gate is imple- 

mented as a single (rectangular-shaped) ;r/2 pulse of duration 
T;j/2 ~ 21 /is with appropriate phase. For calibration simplic- 
ity, a Pauli gate is implemented for ffp G {±0,, ±0,,} as two 
successive 7r/2 pulses, and an identity gate ±7 is implemented 
as an interval of the same duration without the application of 
microwaves. A gate e^^'^°~l^ is implemented as an identity 
gate, but the logical frame of the qubit and subsequent pulses 
are adjusted to realize the relevant change in phase. All pulses 
are separated by a delay of 0.72 /xs. 

To detect the final qubit state, n pulses implement the trans- 
fer 1 1) -> |1,-1) |2,-2) (blue lines in Fig. [T] inset). Two 
additional pulses implement the transfer 1 1) ^ | i) |1, 1) 
(black and red lines in Fig. [T] inset). The ion is then illumi- 
nated for 400 /is by 313-nm light resonant with the cycling 
transition, and the resulting fluorescence is detected with a 
photomultiplier. The entire sequence experiment (from ini- 
tialization through detection) is repeated 100 times (for each 
sequence) to reduce statistical uncertainty. On average, ap- 
proximately 1 3 photons are collected from an ion in the bright 
|2,— 2) state, but only 0.14 are collected from an ion in the 
dark |1,1) state (due largely to laser light scattered from 
the trap surface). To normalize the detection and to elimi- 
nate errors due to slow fluctuations in laser power, each se- 
quence experiment is immediately followed by two reference 
experiments, where the ion is prepared in the | |) and | X) 
states, respectively, and the above detection protocol is imple- 
mented. From the resulting blight and dark histograms [inset 
to Fig.|2|b)], we take the median to establish a threshold for 
1 1) and 1 1) detection. 

Results are shown in Fig. |2] Sequence length refers to the 
number of computational gates in a sequence. We implement 
sequences of lengths 1, 3, 8, 21, 55, 144, 233, 377, 610, and 
987, with 100 different sequences at each length, for a total 
of 1000 unique sequences. With the 21-/is 7r/2 duration used 
here, a sequence of 987 computational gates requires approx- 
imately 64 ms to complete. Our current software limits the 
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FIG. 2: (Color online) Results of the single-qubit benchmarking experiments, (a) Histogram of sequences of a given length with a given 
fidelity. Fidelity is discretized to 0.01 precision because 100 experiments were performed for each sequence, (b) Mean fidelity for each 
sequence length with error bars. The black trace is a least-squares fit to Eq. |lj yielding an EPG of 2.0(2) x 10^^. (Inset) Summed histogram 
of bright and dark calibration experiments with a red line indicating the detection threshold. 



experiment to sequences of length < 1 , 300 gates. 

Theoretically, the average probability for obtaining a cor- 
rect measurement result (the fidelity) after a sequence of 
length I is irOl 



+ -(l-^/if)(l-24)' 



(1) 



where dif describes errors in initialization and measurement 
and S'g is the EPG. A least-squares fit of the observed de- 
cay in fidelity to Eq. ^ yields ^ = 2.0(2) x 10"^ and dif = 
2.7(1) X 10-^. Here, dif is limited by imperfect laser polar- 
ization caused by inhomogeneities in the birefringence of the 
cryogenic windows of the vacuum enclosure. 

The following systematic effects may contribute to the 
EPG: magnetic -field fluctuations, microwave phase and fre- 
quency instability and resolution limits, ac Zeeman shifts, 
pulse amplitude and duration fluctuations, microwave-ion- 
coupling strength fluctuations, decoherence caused by unin- 
tended laser illumination of the ion, and off-resonant excita- 
tion to other levels in the ground-state hyperfine manifold. 

During the benchmarking, we calibrate the qubit transition 
frequency approximately every 60 s. The difference between 
each frequency recalibration and the first calibration is plotted 
in Fig. for the time period corresponding to the data in 
Fig. [2] Monte Carlo simulations of the sequences indicate an 
EPG contribution of 4 = jSA^, where j3 = 1.91 x IO^VHz^ 
and A is the detuning of the microwave frequency from the 
qubit frequency (assumed constant for all of the sequences). 
In the absence of recalibrations, the root-mean-square (rms) 
difference of 25 Hz would give a predicted EPG of 1 .2 x 10^^. 
However, with regular recalibration, the rms difference in fre- 
quency between adjacent calibration points (15 Hz) gives a 
predicted contribution to the EPG of 0.4 x 10^^. The mi- 
crowave frequency and phase resolution are 0.37 Hz and 
1 .5 mrad, respectively, leading to a predicted EPG contribu- 
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FIG. 3: Changes in (a) qubit transition frequency and (b) 7t/2 dura- 
tion during the benchmarking experiments. Change is defined as the 
difference between the recalibrated value and the first calibration. 
Typical transition frequencies and 7t/2 durations are approximately 
1.250 7385 GHz and 20.50 /xs, respectively. 



tion of less than 10 ^. 

A theoretical estimate for the expected ac Zeeman shift of 
the clock (qubit) transition yields a value of less than 1 Hz. 
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In principle, this shift can be determined by comparing the 
qubit frequency measured in a Ramsey experiment with that 
of a Rabi experiment. Such back-to-back comparisons yielded 
values ranging from H-14 Hz to -10 Hz, each with errors of 
approximately 2 Hz. The source of this variation and of the 
discrepancy with theory is not known, but if we assume, as 
a worst-case scenario, a miscalibration of 15 Hz for the fre- 
quency, we estimate an EPG contribution of 0.4 x 10^^. 

One measure of errors caused by qubit-frequency fluctua- 
tions (e.g., from fluctuating magnetic fields) is to character- 
ize decoherence as an exponential decay through a T2 process 
QUOl. To check this, we implement a Ramsey experiment. 
The ion is prepared in \ l), and we apply a single n/2 pulse. 
After waiting for an interval t/2, we apply a n pulse to refocus 
the qubit, and following another interval t/2, we apply an ad- 
ditional n/2 pulse, ideally restoring the qubit back to 1 1). An 
exponential fit of the resulting decay in the 1 1) state probabil- 
ity over periods T < 100 ms gives T2 — 0.38(4) s. Assuming 
this value of T2 also describes frequency fluctuations at times 
on the order of the gate pulses, we predict an EPG contribu- 
tion of 9 X 10^^. Because this exceeds the benchmark value, 
we believe that the noise at shorter periods, in this experiment, 
is smaller than that predicted by a simple exponential fitted at 
longer durations. 

We recalibrate the 7t/2 duration approximately every 120 s 
with a sequence of 256 in-phase n/2 pulses [Fig. [3jb)]. 
Monte Carlo simulations indicate an EPG contribution of 
#g = 7(At)^, where y =2.7 x 10"^/^s^ and At represents 
a miscalibration in the n/2 time (assumed constant for all se- 
quences). In the absence of recalibration, the 23-ns rms drift 
would correspond to anEPGofO.l x 10^^;fromthe estimated 
residual miscalibration between points of 5 ns, we predict an 
EPG contribution of less than 10^^. 

We characterize pulse-to-pulse microwave power fluctua- 
tions by turning on the microwaves continuously and sam- 
pling the power every 10 ns. The integral of the sampled 
power over a 25 -/zs interval is proportional to the total ro- 
tation angle during a pulse of the same duration. We perform 
this integral 12 times, with each 25 -/is interval following the 
previous one by 10 s. Within 120 s after turning on the mi- 
crowaves, we observe a 1% drift in the power. If the pulse- 
to-pulse variation in microwave power is, in fact, this large, 
it corresponds to an EPG contribution of 3 x 10^^. How- 
ever, after a 20-min warm-up interval, we measure a pulse- 
to-pulse power variation of only 0.1%, corresponding to an 
EPG contribution of 0.03 x 10^^. Because the duty cycle of 
the benchmarking experiment is not constant, with sequences 
of different lengths and calibration experiments interspersed 
throughout, it is difficult to assign a specific EPG contribution 
to this effect. However, we do observe larger EPG at higher 
microwave powers, consistent with temperature effects play- 
ing a role at these higher powers. 

To investigate unintended laser light as a source of deco- 
herence, (e.g., from optical pumping), the ion is prepared in 
I D and is allowed to remain in the dark for varying durations. 
We observe no decay in the 1 1) state probability with an un- 



certainty of 2 X 10^^//xs, corresponding to the absence of gate 
errors during the 65-/is randomized gate interval with an un- 
certainty of 1 X 10^^. Similar results are obtained for an ion 
prepared in 1 ■ 

Microwave-induced transitions from the qubit levels into 
other Zeeman levels within the ground-state hyperfine mani- 
fold can be inferred by observing an asymmetry between se- 
quences ending in | I) and those ending in | f). While the 
|2, —2) state fluoresces with 13 photons detected on average, 
other hyperfine states yield, at most, 1.3 photons during the 
400-/is detection period. Therefore, transitions from the qubit 
manifold to other levels would show up as a loss of fidelity 
for sequences ending in \ l), while they would not affect the 
apparent fidelity of sequences ending in | f). For the bright 
sequences in Fig. [2] the EPG is 2.2(5) x 10^^, while for the 
dark sequences it is 2.0(5) x 10^^. We conclude that qubit 
leakage contributes an EPG of < 0.2(7) x 10^^. Similarly, 
if ion heating contributes to the EPG, it should appear as a 
deviation from exponential decay in the benchmarking data, 
which we do not observe. 

For future work, it seems likely that microwave power 
fluctuations could be controlled passively through a suitable 
choice of amplifiers and switching circuitry or actively via 
feedback. Shorter pulses at higher microwave powers would 
diminish errors associated with fluctuating qubit frequency, 
but errors due to off-resonant transitions become more of a 
concern in this regime. Off-resonant transitions could be sup- 
pressed with the use of appropriately shaped pulses, which 
concentrate the microwave spectrum near the qubit transition 
frequency. Self-correcting pulse sequences [21 1 could be used 
to reduce the effects of errors in n/2 duration and transition 
frequency. In a multizone trap array, single-qubit gates imple- 
mented with microwaves will be susceptible to cross talk be- 
tween zones; however, this effect can be mitigated with care- 
ful microwave design, the use of nulling currents in spectator 
zones [16J, and the use of composite pulses L21i . A demon- 
stration of two-qubit gates with errors small enough to enable 
scalable quantum computing remains challenging, but high- 
fidelity single-qubit gates should make this task easier For ex- 
ample, many refocusing and decoupling techniques are based 
on single-qubit gates and can reduce errors during two-qubit 
gate operations L22J . 

This work was supported by lARPA, NSA, DARPA, ONR, 
and the NIST Quantum Information Program. We thank U. 
Warring, M. Biercuk, A. VanDevender, J. Amini, and R. B. 
Blakestad for their help in assembling parts of the experiment, 
and we thank J. Britton, S. Glancy, A. Steane, and C. Bennett 
for comments. This article is a contribution of the U.S. Gov- 
ernment, not subject to U.S. copyright. 



* kbrown@iontrapping.com; present address: Georgia Tech Re- 
search Institute, 400 10th Street Northwest, Atlanta, Georgia 
30318, USA. 



5 



' Present address: QUEST, Leibniz Universitat Hannover, Im 
Welfengarten 1, D-30167 Hannover and PTB, Bundesallee 100, 
D-381 16 Braunschweig, Germany. 

[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, Cam- 
bridge, UK, 2000). 

[2] T. D. Ladd et al.. Nature (London) 464, 45 (2010). 

[3] D. P DiVincenzo, Fortschr. Phys. 48, 771 (2000). 

[4] J. Preskill, Proc. R. Soc. Lond., Ser. A 454, 385 (1998). 

[5] E. Knill, Nature (London) 463, 441 (2010). 

[6] C. A. Ryan, M. Laforest, and R. Laflamme, New J. Phys. 11, 
013034 (2009). 

[7] S. Olmschenk et al.. New J. Phys. 12, 1 13007 (2010). 

[8] J. Benhelm et al., Nat. Phys. 4, 463 (2008). 

[9] C. Ospelkaus et al.. Nature (London) 476, 181 (2011). 
[10] E. Knill et al, Phys. Rev. A 77, 012307 (2008). 



[11] J. F. Poyatos, J. L Cirac, and R ZoUer, Phys. Rev. Lett. 78, 390 

(1997). 

[12] S. Bravyi and A. Kitaev, Phys. Rev. A 71, 022316 (2005). 
[13] E. Knill, Nature (London) 434, 39 (2005). 
[14] M. J. Biercuk et al., Quantum Inf. Comput. 9, 920 (2009). 
[15] J. M. Chow et al., Phys. Rev. A 82, 040305 (2010). 
[16] C. Ospelkaus et al., Phys. Rev. Lett. 101, 090502 (2008). 
[17] L. Deslauriers et al., Phys. Rev. Lett. 97, 103007 (2006). 
[18] J. Labaziewicz et al., Phys. Rev. Lett. 100, 013001 (2008). 
[19] K. R. Brown et al.. Nature (London) 471, 196 (201 1). 
[20] S. Seidelin et al., Phys. Rev. Lett. 96, 253003 (2006). 
[21] M. H. Levitt, Prog. Nucl. Magn. Reson. Spectrosc. 18, 61 
(1986). 

[22] L. Viola, S. Lloyd, and E. Knill, Phys. Rev. Lett. 83, 4888 
(1999). 



